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The Landau theory of phase transitions of Bao.8Sro.2Ti03 thin film under external electric field 
applied in the planar geometry is developed. The interfacial van-der-Waals field E z = 1.1 x 10 8 
V/m oriented normal to the film-substrate interface was introduced in to the model calculation to 
explain experimentally observed behavior of the polarization as a function of planar electric field. 
The E z misfit strain phase diagram of the film is constructed and discussed. 

PACS numbers: 73.61.-r, 77.55.-g, 77.55.Px, 77.55. fe, 77.80. Fm, 77.80.bn 



I. INTRODUCTION 

Thin films of (Bai^Sr^TiOs (BST) are very attrac- 
tive for practical applications in various fields of micro- 
electronics and optoelectronics. Thin films of BST pos- 
sess stable behavior in a wide temperatures range, a high 
dielectric permittivity and small losses that is very impor- 
tant for tunable microwave applications . The struc- 
ture and properties of thin films are well known to exhibit 
a number of deviations from those of bulk ceramics and 
single crystals. The change in the ferroelectric state is as- 
sociated with high internal stresses present in thin films. 
They are due to a substrate-film lattice mismatch, con- 
siderable difference between the linear thermal expansion 
coefficients of the film and substrate, and spontaneous 
strain which appear at the paraelectric-ferroelectric tran- 
sition when the heterostructure is cooled from the deposi- 
tion temperature. A better understanding of correlations 
between various strain factors, which determine proper- 
ties of thin films, is obviously required to achieve reli- 
ability of ferroelectric properties desirable for practical 
applications. 

Interdigital electrodes and parallel-plate electrodes are 
two different types of electrodes used to measure dielec- 
tric properties at different frequencies. At present, many 
studies are focused on the influence of the film/electrode 
interface on the dielectric properties of thin films. In re- 
cent years, great attention was paid to interfacial "dead 
layer" with low permittivity between the film and elec- 
trode 0-0| • These interfacial nano- layers are negligible in 
the case of bulk samples but may be comparable with the 
film thickness. Therefore, interfacial "dead layers" may 
seriously influence the effective dielectric response per- 
pendicular to the thin film plane when usual plane capac- 
itor geometry metal-dielectric-metal (MDM) is used for 
dielectric measurements with a small distance between 
top and bottom electrodes. 

It is worth mention that interfacial "dead layers" are 
of minor influence for the in-plane direction of the ex- 
ternal field created by planar electrodes deposited onto 
film surface where the gap between the electrodes is not 



connected with the film thickness. The influence of the 
" dead layer" is maximized in the MDM case and negligi- 
ble when planar electrodes are used. However in the lat- 
ter case one-dimensional phenomenological models may 
yield wrong results. Models describing a change of the 
polarization direction under external field similar to that 
developed for the spin-flop transitions in ferromagnetic 
materials Q are obviously required. 



In this paper we developed the phenomenological 
model for epitaxial Bao.sSro^TiOa (BST08) thin film on 
cubic substrate with the aim to describe polarization be- 
havior under external electric field applied along the film 
plane. Utilizing the results of previous theoretical work 
, we construct a misfit strain ( u m )- electric field ( E ) 
phase diagram BST08 assuming presence of two factors: 
misfit strain imposed by the substrate and electric field 
induced by the film-substrate interface. Existence of this 
interfacial field follows form the absence of the symmetry 
plane at the film-substrate interface. The physical origin 
of interfacial field is discussed in the last section of the 
paper. Diffuse character of phase transitions in ferroelec- 
tric thin films (dielectric maxima is always broad) is an 
indirect evidence of the presence of such field. It is worth 
note that phase transitions in ferroelectric materials un- 
der external filed conjugated with the order parameter 
exhibit quite similar diffused character. 



In the frame of the proposed model, we calculated the 
dependence of each component of the polarization vec- 
tor as function of planar field applied parallel to the film 
surface. Experimental measurements of polarization hys- 
teresis loops were performed for the 40-nm-thick BST08 
film deposited onto (001) MgO substrate. Using the the- 
oretical model and the lattice parameters of this film we 
calculated the polarization as a function of planar field 
and estimated the value of the interfacial field. 
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II. PHENOMENOLOGICAL THEORY OF 
POLAR STATES IN FERROELECTRIC THIN 
FILM UNDER EXTERNAL ELECTRIC FIELD 

Following Pertsev et al [l(| the thermodynamic phe- 
nomenological potential of a thin film is constructed us- 
ing the thermodynamic potential for the bulk sample. 
The thermodynamic phenomenological potential BST 
solid solutions and phase diagrams for BST films epitax- 
ially grown onto (001) surface of a cubic substrate were 
recently reported Q. According to this work, only one 
ferroelectric order parameter -polarization can be consid- 
ered for BST08 thin film. The relevant expression for the 
cubic perovskite film deposited onto (001) surface of the 
cubic substrate is given as 

G = ai(pl+p 2 y ) + a 3 p 2 z + an(px + P%) + a 33Pt 
+ai2P 2 x p 2 y + a 13 (p 2 x + p 2 y )p 2 z + am(pl + p 6 y + p 6 z ) 

+ ail2[pi{p 2 y + Pi) + Py(pl + Pi) + PUPI +P 2 y)\ 

+ai23P 2 x pl,p 2 z + a nn(pl + P* y + pI) (1) 
+a 1112 [pl{p 2 y +p 2 z ) +pl(p 2 x + p 2 z ) + p 6 z (p 2 x + p 2 y)\ 

i / 4 4 i 4 4 i 4 4\ 

+ ail22(P X Py + P X P Z + PyP Z ) 

i / 4 2 2 i 242, 22 4\ 

+ ail23{P X PyP Z + P X PyP Z + P X PyP Z )- 

All coefficients of the potential (1) were calculated for 
particular concentration BST08 using the coefficients of 
the end-members of the solid solution according to the 
formalism developed for BST solid solution Q . The high- 
temperature asymptotic values of the coefficients in the 
quadratic terms of the potential for BST08 yield 

at = 4.791 x 10 5 (T - 328) - 0.857 x 10 10 u m Jm/C 2 , 
a 3 = 4.791 x 10 5 (T - 328) + 1.014 x 10 10 u m Jm/C 2 , 

where u m is misfit strain. All calculated coefficients 
are listed in Table I. 

Strain of a film with respect to the bulk BST08 sample 
are 

6 1 U m ; 

62 = Um, 

e 3 = -0.00659p 2 - 0.00659p 2 + 0.04276p 2 

- 0.84668u m , (2) 
e 4 = 0.02673p y p z , 
e s = 0.02673?^, 
e 6 = 0. 

In the following we consider the single-domain film de- 
posited onto (001) face of the cubic substrate with z 
axis being normal to the film surface. External electric 
field E x is applied using planar electrodes deposited on 
the film surface and aligned parallel to one of the sub- 
strate's cubic axes - y. The interfacial van-der-Waals 



TABLE I: Coefficients of the phenomenological potential (1) 
for BST08 film deposited onto (001) surface of the cubic sub- 
strate 



Coefficient 




BblOS 


Units 




4.791(T 


- 328) - 0.857 x 10 5 u m 


xlO 5 Jm/C 2 


«3 


4.791(T 


- 328) + 1.014 x 10 5 it m 




(111 




1.6837 


xl0 8 Jm 5 /C7 4 


(ll> 




4.1796 




ai3 




4.4478 




ffl33 




-1.4913 




ttlll 




1.0352 


xl0 9 Jm 9 /C7 6 


dll2 




-1.5600 




ffll23 




-0.6080 




ami 




3.0904 


x 10 10 Jm 13 /C 8 


ani2 




2.0232 




CH122 




1.3096 




aii23 




1.0936 





field E z arises at the interface of two surfaces and de- 
pends strongly on the technological parameters of the 
film growth. To analyze the effect of the external electric 
field, additional terms —E x p x — E z p z should be included 
into potential (1). One more term responsible for the de- 
polarizing filed should be included in the potential (1). 
However, in our simple model we assume that inhomoge- 
neous polarization on the film surface equal to zero due 
to presence of inhomogeneous interfacial van-der-Waals 
field. 

Fig. 1 shows the misfit strain ( u m ) - interfacial 
field E z phase diagram under zero external planar field 
( E x = ). A component of polarization p z normal to 
the film surface exists at any u m . Two additional field- 
induced phases appear at positive u m values with addi- 
tional component of the polarization along the cubic axis 
(Px p z ) and along the diagonal of the film surface (p x 
Px p z )- The phase transition between (p x p z ) and p x p x 
p z ) is of first-order type (bold line in Fig.l), while tran- 
sition between (0 p z ) and (p x p z ) is of second-order 
type (dashed line in Fig.l). At E z = 0, the phase transi- 
tion between (p x p x p z ) and (0 p z ) (first order) occurs 
at u m — 0.23 x 10~ 3 , however, (p x p x p z ) is metastable 
up to u m = -0.05 x 10~ 3 . 

Dotted lines in Fig. la show a range of 180° domains 
of the tetragonal film while above and below these lines 
the film is in single-domain state with alternative direc- 
tion of polarization +p z and — p z , respectively. Enlarged 
sections of the diagram near E z — and u m — is shown 
in Fig. lb. According to this diagram, on increasing u m 
at E z — the polarization changes from the direction 
orthogonal to the substrate (0 p z ) to the tilted one - 
{Px Px p z )- Note that polarization direction (p x p z ) is 
possible at nonzero E z values. 

In the following we analyze the polarization of the film 
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FIG. 1: a) Misfit strain u m - interfacial field E z phase diagram BST08 thin film, b) Enlarged sections of the diagram near 
E z = and u m — 0. Solid and dashed lines correspond to first- and second-order phase transitions, respectively. Thin dotted 
lines correspond to the metastable states. 



as a function of external electric field E x . In various parts 
of the phase diagram dependence of polarization is quite 
different especially near the phase transition lines, and 
very complicated around the multiphase point. 

Polarization as a function of ac external electric field 
and interfacial field E z = const was calculated using dy- 
namic equations of Landau-Khalatnikov with a uniform 
relaxation time for different component of the polariza- 
tion. 



_ dp x 
dt 
dpy 
dt 

pVz 
dt 



dG 
dp x 
dG 

dG 
dp z ' 



(3) 



A Runge-Kutte method was employed for numerical 
calculations. Fig. 2 shows typical loops calculated for all 
components of polarization at E z — 10 8 V/m in several 
point of the phase diagram corresponding to the differ- 
ent values of misfit strain. The chosen value of E z used 
in our calculations will be justified below. Dependences 
of all components of polarization as a function of exter- 
nal electric field E x become complicated with increasing 
u m - In the tetragonal phase (0 p z ) the p x component 
(solid line in Fig. 2a) exhibits nonlinear behavior. On 
approaching the phase transitions lines the slop of the 
p x component increases and hysteresis starts to appear 
(Fig. 2b). Finally, in the phase (p x , p x , p z ) the p y com- 
ponent exhibits dome-like feature around zero E x values 
(Fig. 2c) and hysteresis behavior appears at the edges of 
the dome and near zero E x values as well. It worth not- 
ing that all three components of the polarization exhibits 
hysteresis inside the {p x , p x , p z ) phase. 



III. EXPERIMENTAL RESULTS 

The transparent and mirror-smooth heteroepitaxial 
BST08 thin film (thickness 40 nm) was deposited on 
(OOl)MgO single crystalline 0.5 mm-thick substrate by 
rf sputtering of polycrystalline target of the correspond- 
ing stoichiometric composition. The temperature of the 
substrate was kept at 650 — 750°C. Details of the growth 
conditions have been previously reported [Tl|] . The epi- 
taxial relationships between the film and the MgO sub- 
strate have been confirmed by X-ray diffraction (XRD): 
(001)bst II (001) M gO and [100] BS t || [100] Mg o- XRD 
study revealed that the room-temperature unit cell pa- 
rameter of the film normal to the substrate (c = 0.39936 
nm) is smaller than that in plane of the film (a — 0.40032 
nm) . At 600° the lattice parameters of this film were de- 
termined to be c = 0.40118 nm and a — 0.40287 nm. 

For application of external electric field, a periodic ar- 
ray of 440 interdigital Al electrodes were deposited on 
the top of the film using standard optical photolithog- 
raphy. The electrodes thickness, length and width were 
300, 130 and 2 /zm, respectively. The spacing between 
the electrodes was 1 /im. The electrodes were aligned 
parallel to one of the substrate's cubic axes. An optical 
microscope "Leitz Latimet" was used to verify topology 
and quality of the electrodes. Capacity, conductivity and 
leakage current were measured at 1 kHz using Keihtley 
4200SCS analyzer and PM-5 MicroTec probe station. 

Room-temperature polarization-electric-field hystere- 
sis measurements were performed in a frequency range 
10~ 2 - 10 3 Hz using a modified Sawyer- Tower circuit 1131 . 

To estimate the film strain we used the equation 13]: 



u o - ("/ -&s)(T - T ) (4) 

a f 



4 



P (C/m 2 ) P (C/m 2 ) 




a) E = 

FIG. 2: Three components of polarization as a function of ac 

electric field E x under various values of misfit strain at fixed FIG. 3: Three components of polarization as a function of ac 
interfacial field E x = fO 8 V/m. A solid line - p x , dotted line electric field E x under various values of a interfacial field at 
- p y , a dashed line - p z . Figures a-c corresponds to the points fixed misfit strain u m = 0.16 x fO -3 . A solid line - p x , dotted 
indicated in Fig. fa. line - p y , a dashed line - p z . Figures a-c corresponds to the 

points 1-3 indicated in Fig. fa. 
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p x (C/m 2 ) 




-0.2 - 

FIG. 4: Experimental dependence (open symbols) of p x com- 
ponent of the polarization as a function of ac electric field E x 
and theoretical dependence (bold line) of this component at 
E„ = 1.1 x 10 s V/m and u m = 0.16 x 1CT 3 . 



where a/ , a s -are the thermal expansion coefficients of 
the film and substrate, respectively; a/ and a s - are the 
lattice parameters of the film and substrate at To ; uo - is 
initial strain at at To. At room temperature a s =0.4212 
and thermal expansion coefficient a s = 1.2 x 10~ 5 1/K 
above room temperature up to 600°C [HJ, therefore 
a,,=0.4247nm at 600°C. The lattice parameters of bulk 
BST08 solid solution were calculated according Eqs. (7) 
and (8) from Ref. [9]. Using the structural informa- 
tion about end-members of the solid solution [l^, [T(| , we 
found a f = 0.922 x 1CT 5 1/K and a/=0.3998 nm. As- 
suming that spontaneous polarization is zero at 600°C 
and using the relation c/a= (1 + e 3 )/(l + ei), Eqs. (2) 
and (4) at T = T we found u = 2.28 x 1(T 3 . Therefore, 
the room-temperature misfit strain calculated using Eq. 
(4) is u m = 0.16 x 10~ 3 . 

Fig. 3 shows typical loops calculated for all compo- 
nents of polarization at fixed misfit strain u m = 0.16 in 
several points of the phase diagram corresponding to the 
different values of the interfacial field E z . These three 
points are indicated by numbers 1-3 in Fig. la. 

The polarization as a function of ac electric field E x 
calculated at zero interfacial field (Fig. 3a ) does not 
agree with the experimental loop presented in Fig. 4. 
Reasonable agreement can be achieved at E z = 1.1 x 10 8 
V/m. The corresponding theoretical dependence of p x 
component as a function of ac electric field E x is shown 
in Fig. 4 (bold line). 



IV. DISCUSSION AND CONCLUSIONS 

A phenomenological model including three compo- 
nents of the polarization was developed to explain the 



dependence of the in-plane component of the polariza- 
tion p x as a function of ac electric field E x applied paral- 
lel to the film plane. The misfit strain u m of the film was 
found to be equal to 0.16 x 10~ 3 . The interfacial van- 
der-Waals field E z = 1.1 x 10 s V/m oriented normal to 
the film-substrate interface was introduced in the model 
calculation to explain experimentally observed behavior 
of the polarization as a function of planar electric field. 
In the following we justify the existence of the interfacial 
van-der-Waals field. 

A short-range van-der-Waals interaction always exits 
between two atoms due to mutual polarization of elec- 
tronic shells. In some cases, the interaction energy of 
two planar surfaces with the distance r is ~ r~ 2 for unit 
area of one surface interacting with an infinite area of an- 
other surface [l7| . As known from the tip-surface curves 
in atomic-force microscopy pjl the radius of the inter- 
action varies from several angstroms up to several tens 
of nanometers. These van-der-Waals fields are only par- 
tially compensated when two planar surfaces are in con- 
tact. The resulting energy and direction of this filed at 
the interface of two surfaces is a subject of quantum- 
mechanical calculations. In fact, the interfacial van-der- 
Waals field is a result of the polarization induced at the 
substrate. 

If polarization at the external film surface is not equal 
to zero, depolarizing field appears to reduce the spon- 
taneous polarization The direction of the depolar- 
izing field is opposite to the direction of the interfacial 
van-der-Waals field, therefore the depolarizing field re- 
duces the effect of the interfacial one. The interfacial 
van-der-Waals field is a long-range inhomogeneous field 
[l7j , therefore our homogeneous phenomenological model 
should be considered as approximation of the averaged 
inhomogeneous model in which boundary conditions on 
the external surface of the film allows absence of the de- 
polarizing field. 

In general case, the interfacial van-der-Waals field fluc- 
tuates and its direction is not defined. However, under 
external field the van-der-Waals field acquires particu- 
lar direction since fluctuating dipoles become oriented 
under external field. Therefore, below the ferroelectric 
phase transition temperature, the interfacial van-der- 
Waals field has preferred direction normal to the sub- 
strate. 

In fact, the introduction of the interfacial van-der- 
Waals field is equivalent to the introduction of so-called 
"passive layer" [20|. Very likely, that origin of the "pas- 
sive layer" is the effect of the van-der-Waals fields in these 
(the interfacial) regions. The size of the "passive layer" 
depends on the physical properties of the contacting ma- 
terials. Namely, the field strength and its long-range 
properties depend mainly on the spectral characteris- 
tics of the dielectric constants of the contacting materials 

The calculated p x (E x ) dependence (bold line in Fig. 
4) implies that experimentally observed losses (loops in 
Fig. 4) appear due to relaxations in the regions close to 
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electrodes and not due to the losses inside the material of 
the film studied. Therefore, in functional devices based 
on this film one can expect low losses in microwave fre- 
quency range. As follows from Fig. 2, by increasing the 
misfit strain one can increase the slope of the hystere- 
sis loop and therefore increase the tunability. However, 
large positive misfit strains induce the hysteresis behav- 
ior implying large losses inside the material of the film. 
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